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An attitude and orbit control system (AOCS) is developed for a 160-m, 450-kg solar sail 
spacecraft of the Solar Polar Imager (SPI) mission. The SPI mission is one of several Sun- 
Earth Connections solar sail roadmap missions currently envisioned by NASA. A reference 
SPI sailcraft consists of a 160-m, 150-kg square solar sail, a 250-kg spacecraft bus, and 
50-kg science payloads, The 160-m reference sailcraft has a nominal solar thrust force of 
160 mN (at 1 AU), an uncertain center-of-mass/center-of-pressure offset of ±0.4 m, and a 
characteristic acceleration of 0.35 mm/s 1 2 . The solar sail is to be deployed after being placed 
into an earth escaping orbit by a conventional launch vehicle such as a Delta II. The SPI 
sailcraft first spirals inwards from 1 AU to a heliocentric circular orbit at 0.48 AU, followed 
by a cranking orbit phase to achieve a science mission orbit at a 75-deg inclination, over 
a total sailing time of 6.6 yr. The solar sail will be jettisoned after achieving the science 
mission orbit. This paper focuses on the solar sailing phase of the SPI mission, with 
emphasis on the design of a reference AOCS consisting of a propellantless primary ACS and 
a microthruster-based secondary (optional) ACS. The primary ACS employs trim control 
masses running along mast lanyards for pitch/yaw control together with roll stabilizer bars 
at the mast tips for quadrant tilt (roll) control. The robustness and effectiveness of such 
a propellantless primary ACS would be enhanced by the secondary ACS which employs 
tip-mounted, lightweight pulsed plasma thrusters (PPTs). The microPPT-based ACS is 
mainly intended for attitude recovery maneuvers from off-nominal conditions. A relatively 
fast, 70-deg pitch reorientation within 3 hrs every half orbit during the orbit cranking 
phase is shown to be feasible, with the primary ACS, for possible solar observations even 
during the 5-yr cranking orbit phase. 

I. Introduction 

Solar sails are envisioned as a propellantless, high-energy propulsion system for future space exploration 
missions. NASA’s future missions enabled by solar sail propulsion include the Solar Polar Imager (SPI), 
Ll-Diamond, Particle Acceleration Solar Orbiter (PASO), and Interstellar Probe, which are the Sun-Earth 
Connections (SEC) solar sail roadmap missions. 1-3 In particular, the SPI mission is currently being further 
studied by NASA/JPL, and it is now called the SPI Vision mission. 4 Our current understanding of the Sun 
is limited by a lack of observations of its polar regions. The SPI Vision mission utilizes a large solar sail to 
place a spacecraft in a 0.48-AU heliocentric circular orbit with an inclination of 75 deg. Viewing of the polar 
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Figure 1. A solar sail mast with a trim ballast mass (running along a lanyard tape) and tip-mounted roll 
stabilizer bars attached to sail panels. 


regions of the Sun provides a unique opportunity to more fully investigate the structure and dynamics of its 
interior, the generation of solar magnetic fields, the origin of the solar cycle, the causes of solar activity, and 
the structure and dynamics of the corona. 

The SPI mission consists of the initial cruise phase (1.6 yr) to a 0.48- AU circular orbit, the cranking 
orbit phase (5 yr), and the science mission phase (2 yr). A 160-m, 450-kg solar sail spacecraft is considered 
for such a solar sailing mission. 4 A Delta II launch vehicle is able to inject the 450-kg SPI spacecraft into 
an earth escaping orbit with C 3 = 0.25 km 2 /s 2 , and then the sail is to be deployed. The SPI sailcraft first 
spirals inwards from 1 ATJ to a heliocentric circular orbit at 0.48 AU, then the cranking orbit phase begins 
to achieve a 75-deg inclination. The solar sail will be jettisoned after achieving the science mission orbit, 
and the total sailing time is approximately 6.6 yr. 

This paper will focus on the attitude and orbit control of a large sailcraft during its solar sailing phase. A 
reference SPI sailcraft selected for AOCS design consists of a 160-m, 150-kg solar sail, a 250-kg spacecraft bus, 
and 50-kg science payloads. This 160-m reference sailcraft has: a nominal solar thrust force of 160 mN (at 
1 AU), an uncertain center-of-mass/center-of-pressure offset of ±0.4 m, and a characteristic acceleration of 
0.35 mm/s 2 . A baseline AOCS architecture proposed for the SPI sailcraft consists of a propellantless primary 
ACS and a secondary ACS. The primary ACS employs trim control masses running along mast lanyards 
for pitch/yaw control together with roll stabilizer bars at the mast tips for quadrant tilt (roll) control. The 
secondary ACS employs tip-mounted, lightweight pulsed plasma thrusters (PPTs). MicroPPT modules are 
mounted at the tips of the sail booms providing a large control moment arm (112 m). The microPPT-based 
secondary ACS is mainly intended for attitude recovery maneuvers from off-nominal conditions. For more 
general sailing missions, where the sail is never jettisoned, this secondary ACS provides a lower-cost, lower- 
mass propulsion for deployment control and greater redundancy than any traditional reaction-jet control 
system. 

Either AOCS is able to perform a 70-deg pitch reorientation within 3 hrs every half orbit during the 
orbit cranking phase. Such a relatively fast, orbit cranking maneuver is currently being considered by the 
SPI Vision mission study team for possible solar observations even during the 5-yr cranking orbit phase. 

The remainder of this paper is outlined as follows. Section II describes the baseline AOCS architecture 
proposed for the 160-m SPI sailcraft. Section III describes an overview of the SPI Vision mission during its 
solar sailing phase. Section IV describes the attitude-trajectory coupled equations of motion of a sailcraft in 
a heliocentric orbit for attitude/orbit control design. Section V presents trajectory optimization results for 
the SPI Vision mission. Section VI presents the propellantless primary ACS design. 
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Figure 2. An integrated attitude and orbit control system (AOCS) architecture proposed for the 160-m, 450-kg 
SPI sailcraft. 


II. An Integrated AOCS Architecture of the SPI Sailcraft 

A baseline AOCS, consisting of a propellantless primary ACS and a microPPT-based secondary ACS, is 
proposed for a 160-m reference sailcraft of the SPI Vision mission. The primary ACS employs trim control 
masses (TCMs) running along mast lanyards for pitch/yaw control together with roll stabilizer bars (RSBs) 
at the mast tips for quadrant tilt control, as illustrated in Fig. 1. The robustness of such a propellantless 
primary ACS would be further enhanced by a secondary ACS utilizing tip-mounted, lightweight PPTs, which 
are also illustrated in Fig. 1. Such an optional microPPT-based ACS can be employed for attitude recovery 
maneuvers from off-nominal conditions as well as for a spin-stabilized safe mode. It can also be employed as a 
backup to the conventional ACS of the spacecraft prior to sail deployment and during pre-flight sail checkout 
operation, if necessary. A conventional bus ACS is required for the SPI mission as the sail is jettisoned at 
the start of the science mission phase. As an alternative to conventional approaches, the microPPT-based 
ACS option promises lower mass, lower cost, and greater redundancy. 

The attitude determination system (ADS) is a critical subsystem of most spacecraft AOCS. An ADS 
of particular interest for solar sail applications is the Inertial Stellar Compass (ISC) recently developed 
by Draper Laboratory for an NMP ST6 flight validation experiment. 6 The ISC is a miniature, low-power 
ADS developed for use with iow-cost microsatellites. It is suitable for a wide range of future solar sail 
missions because of its low-mass, low-power, and low-volume design and its self-initializing, autonomous 
operational capability. The ISC is composed of a wide field-of-view active-pixel star camera and microgyros, 
with associated data processing and power electronics. It has a total mass of 2.5 kg, a power requirement of 
3.5 W, and an accuracy of 0.1 deg (lcr). It is planned to be flight validated within few years. Some recent 
advances in microsatellite technologies, including the ISC, need to be exploited to complete an integrated 
low-cost, low-risk, low-mass, low-power, and iow-voiume AOCS for sailcraft. 

The proposed AOCS architecture with various ADCS options, including a conventional ADCS for the 
spacecraft, is illustrated in Fig. 2. Detailed analysis and design of this baseline AOCS architecture, as 
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applied to a fight validation mission of a 40-m solar sail in a dawn-dusk sunsynchronous orbit, can be found 
in Refs. 5,7-9. Other ACS options, employing a control boom or control vanes, can be be found in Ref. 10. 
A potential control-structure interaction caused by the significant torsion and bending of flexible masts of 
a large sailcraft is discussed in Ref. 11. The SPI mission architecture and its 160-m sailcraft equipped 
with the proposed baseline AOCS will be applicable with minimal modifications to a wide range of future 
solar sail flight missions with varying requirements and mission complexity, including a solar sailing mission 
for intercepting, impacting, and deflecting near-Earth asteroids . 12,13 More details of the microPPT-based 
secondary ACS proposed for the SPI sailcraft can be found in Ref. 14. 

III. SPI Solar Sailing Mission Overview 

The SPI mission basically consists of the initial cruise phase (1.6 yr) to a 0.48-AU circular orbit, the 
cranking orbit phase (5 yr), and the science mission phase (2 yr). The initial cruise phase of spiraling 
toward the sun using a simple in-plane steering law with a fixed 35-deg sun angle is illustrated in Fig. 3 
for a reference sailcraft with a characteristic acceleration of a c — 0.35 mm/s 2 (corresponding to a lightness 
number of A = 0.06). The characteristic launch energy, which is defined as the square of the hyperbolic 
excess speed v^, was assumed as C 3 = 0.25 km 2 /s 2 for a Delta II launch vehicle. The initial position was 
simply assumed as (X,Y) = (1,0) AU. In Fig. 3, the sailcraft speed is compared to the circular orbit speed 
(dotted red line) at the corresponding orbital position r of the sailcraft. Although the use of this simple 
steering law does not exactly result in a circular orbit at 0.48 AU, it is presented here as a simple reference 
trajectory. Various optimization-based trajectories of the cruise phase will be discussed later in this paper. 

The cranking orbit phase for achieving a 75-deg inclination by employing a simple out-of-plane steering 
law of ±35-deg sun angle change every half orbit is illustrated in Fig. 4. The elliptic orbital nature of the 
cranking orbit phase is due to the radial component of the solar radiation pressure force. 

For typical solar sailing flight optimization problems, trajectory models which are decoupled from attitude 
dynamics were often used in the past . 15-17 However, the effect of attitude motion of large solar sails on the 
solar sailing trajectory is of current practical concern. Consequently, the development of a six-degree-of- 
freedom, orbit-attitude coupled dynamical model of solar sail spacecraft for trajectory optimization and 
simulation is of current practical importance . 18,19 

In the next section, a simple attitude-trajectory coupled dynamical model of solar sail spacecraft is 
described for attitude control design and integrated orbit-attitude control simulation. A quaternion-based 
attitude control system approach, previously developed in Refs. 20-22, is proposed for the AOCS design of 
the SPI sailcraft because of its simplicity of accommodating the thrust vector control command. 

IV. Coupled Attitude and Orbital Dynamics 

The orientation of the solar sail thrust vector is often described in terms of the cone and clock angles . 15-19 
These two angles are the typical control inputs used in solar sailing trajectory optimization. However, there 
are at least two different sets of cone/clock angles selected for trajectory optimization in the literature. 
In this section, the trajectory-attitude coupled equations of motion, employing such two different sets of 
cone/clock angles, are briefly described for the purposes of trajectory design, attitude control design, and in- 
tegrated orbit-attitude control simulation. Detailed discussions of such trajectory-attitude coupled equations 
of motion can be found in Ref. 19. 

A. Cone and Clock Angles 

Let {/, J, K } and {r, U, <f>} be respectively a set of right-handed, orthonormal vectors of the heliocentric 
ecliptic rectangular and spherical coordinate reference frames, as illustrated in Fig. 5. The two angles, ijj 
and <j>, are called the ecliptic longitude and latitude of the sailcraft position, respectively; 0 < ip < 360 deg 
and —90 deg < <j> < +90 deg. 

The sailcraft position vector is then expressed as 

r — rf 

= (r cos 4> cosi/Of + (rcos<f>sini p)j + (rsin<p)K 

= XI + YJ + ZK (1) 
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Figure 4. The cranking orbit phase for achieving a 75-deg inclination. 


where r = |f] is the distance from the sun to the sailcraft. 

The orientation of a unit vector normal to the sail plane, h, is described in terms of the cone angle a and 
the clock angle /?, illustrated in Fig. 6, as follows: 


n = (cos a)f + (sin a sin /3 )ip + (sinacos/3)0 


(2) 


where 


cos/3 = 


cos a = r ■ n 
f x (n x f) 


■ <i> 


|f x (n x f)| 

0 < a < 90 deg 
0 < j8 < 360 deg 

As also illustrated in Fig. 6, the sailcraft body-fixed basis vectors {1,2,3} are assumed to be aligned with 
when a = 0 = 0, and the sailcraft roll axis is defined to be perpendicular to the sail surface; i.e., 

i = n. 

A different set of the cone and clock angles, (a, S), shown in Fig. 7, can also be employed as follows. 
Let {r,6,k} be a set of basis vectors of an osculating orbital plane, as illustrated in Fig. 8. A body-fixed 
rotational sequence to {1,2,3} from {/, J,K] is then described by successive coordinate transformations of 
the form 21 
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Figure 5. Heliocentric ecliptic rectangular coordinates ( X,Y,Z ) and spherical coordinates (r, ip, rp). 


The orientation of a unit vector normal to the sail plane, h , is then described in terms of a and <5, as 
illustrated in Fig. 7, as follows: 


n = (cosa)r + (sin a sin 5)0 + (sinacos5)fc 


and 


cos 5 = 


cos a = r ■ n 
f x (n x f) 


| f x (fix r)| 
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0 < 8 < 360 deg 

The two different clock angles, 0 and 5, are related as 
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B. Solar Radiation Pressure Force 


An ideal model of the solar radiation pressure (SRP) force is considered here. The SRP force vector (per 
unit mass) acting on the sailcraft is then described in various coordinates as follows: 


where 


F = Fo(^ • fifh 


= F r f + F^ip + F^ip 


= Rf + TO + Nk 


= Fxl + Fyj + F Z K 

(6) 

F° = { T i)C< 

(7) 
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Figure 6. Cone angle a, clock angle 0, and sailcraft orientation when a = 0 = 0. 



Figure 7. Cone angle a and clock angle 6 (Ref. 15). 


where r® = 1 AU = 149597870.691 km is the mean distance from the sun to the earth, and a c is the 
characteristic acceleration of the sailcraft at 1 AU defined as 


&c — 


Ei 

m 


2r)PA 

m 


= 77(9.1 x W~ 6 )A/m 


(8) 


and F s is the nominal solar thrust force at 1 AU for an ideal flat solar sail with a reflectivity coefficient 77 
and a sail area of A. 

Furthermore, we have 


F r 
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Figure 8. Orbital geometry (illustrated for a near-circular orbit). 


C. Gauss’s Form of the Variational Equations (Osculating Orbital Elements) 

Various forms of trajectory equations are available, as described in Ref. 19. For the purpose of preliminary 
trajectory design, attitude control design, and integrated orbit-attitude control simulation, Gauss’s form of 
the variational equations in terms of osculating orbital elements is chosen here although this set of equations 
have a singularity problem when i — 0. 

A set of six first-order differential equations, called Gauss’s form of the variational equations, in terms of 
osculating orbital elements, are given by 23 


where 
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D. Attitude Dynamics and Control Model Employing Quaternions 

A quaternion-feedback attitude control system approach, previously developed in Refs. 20-22, is proposed for 
solar sail attitude applications because of its simplicity of accommodating the commanded cone and clock 
angles of the solar thrust vector into a solar sail attitude control system. 
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Euler’s attitude dynamical equations of motion of a rigid sailcraft axe given by 


I\Co\ - (I2 - I 3 )u>2u 3 = u\ + di ( 12 a) 

I2UJ2 — (I 3 ~ A) w 3 Wl = U2 + d,2 ( 12 b) 

Is<^3 — (A — l2)uiU>2 = U 3 + d 3 ( 12 c) 

where (aq , u > 2 , w 3 ) are the angular velocity components, (h , I 2 , I 3 ) the principal moments of inertia, (■ ui , u 2 , u 3 ) 
the attitude control torques, and ( di,d.2,d 3 ) the disturbance torques. 

The kinematic differential equations in terms of attitude quaternions are described as 
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where (qi, 92, 93, 94) are the inertial attitude quaternions which are the typical outputs from an attitude 
determination system as illustrated in Fig. 2. 

The quaternion-feedback attitude control logic, proposed for solar sailing applications, is simply a PID 
control logic of the form 19-21 
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where (ei,e2,e 3 ) are the roll, pitch, and yaw components of attitude-error quaternions (e 1 ,e 2 ,e3,e 4 ), and 
(Ki, Tt, Di) are PID control gains to be properly determined. 

A saturation control logic of accommodating the actuator torque and slew rate constraints is also given 
by 


Ui = - sat { Ki sat[e< + — f e.dt] + D { w< 1 ; i = 1 , 2, 3 
O m .x { Li T t J J 


( 15 ) 


and the variable limiter L t is self-adjusted as 


Li ^ min j \Z2u, | e j | , cu max ^ 


( 16 ) 


where w max is the maximum slew rate constraint for each axis (if required) and a t is the maximum angular 
control acceleration. Details of this nonlinear PID control logic which is applicable to solar sail attitude 
control can be found in Refs. 21 - 22 . 

The attitude-error quaternions are computed using the desired or commanded attitude quaternions 
(?1C, 92c, 93c, 94c) and the actual inertial attitude quaternions (91,92,93,94), as follows: 
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The desired attitude quaternions (910 92c, 93c, 94c), corresponding to the desired solar thrust vector di- 
rection described by (a, /3), are given by 
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cos(0/2) 0 sin(0/2) 0 
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where 0 and <j> axe the ecliptic longitude and latitude of a sailcraft inertial position, respectively, as shown 
in Fig. 5. 

The desired attitude quaternions ( <7ic,92c,<?3c>94c ). corresponding to the desired solar thrust vector di- 
rection described by (a,<5), are given by 
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where i, w, f2, and 6 are the various orbital angles shown in Fig. 8. 


V. Solar Sailing Trajectory Optimization for the SPI Mission 

Some preliminary study results of developing solar sail trajectory optimization tools, as applied to the 
SPI mission, are presented here. Optimization approaches under study are indirect- method global optimizers 
such as downhill simplex, genetic algorithms, and simulated annealing. 

Figure 9 shows a baseline, optimal trajectory design by Carl Sauer at JPL for achieving a circular orbit 
at 0.48 AU with a 75-deg inclination. The monotonically decreasing semimajor axis and the corresponding 
variation of eccentricity can be seen in this figure during the initial cruise phase to 0.48 AU. The eccentricity 
remains constant during the orbit cranking phase. The eccentricity is finally nulled after cranking is complete. 
The sail angles (a*, 5') in Fig. 9 are the optimal cone and clock angles, respectively. However, the clock 
angle shown in Fig. 9 is actually /3 defined in Fig. 6, not <5 as defined in Fig. 7. Note that the cone angle, a, 
is nearly kept constant at 35 deg throughout the sailing mission. 

Simple sail-steering laws for maximizing the rate of change of the semimajor axis and the inclination, 
described in Ref. 15, were applied to the SPI mission. These sail-steering laws are based on the variational 
orbit equations to optimally change individual orbit elements, and they provide the optimal cone and clock 
angles (a*, 6*) in closed, analytical form. Figure 10 shows the result of applying these simple sail-steering 
laws, given by Eqs. (4.56), (4.63), and (4.70) of Ref. 15, to the SPI mission. The switch from orbit radius 
reduction to the cranking phase occurs once the sail reaches the target semimajor axis of 0.48 AU. The 
simulation used a nominal characteristic acceleration of 0.312 mm/s 2 (at 1 AU), an 8-deg initial inclination, 
and C 3 = 0.25 km 2 /s 2 . The sailcraft achieves the desired semimajor axis and inclination but the final orbit 
is eccentric while the desired SPI mission orbit is circular. The corresponding three-dimensional orbital 
trajectory is illustrated in Fig. 11. Note that the semimajor-axis control law described in Ref. 15 was 
not intended to be used for controlling both the semimajor axis and the eccentricity. Such a simple sail- 
steering law in closed, analytical form was examined as the first step towards developing solar sail trajectory 
optimization tools. 

In an attempt to control the semimajor axis and the eccentricity simultaneously, we have applied the 
indirect methods to the initial cruise phase for achieving a circular orbit at 0.48 AU. Figure 12 shows the result 
of applying the semimajor axis change optimization using Downhill-Simplex, a global optimization method 
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Figure 9. A baseline, optimal trajectory design by Carl Sauer at JPL (C 3 = 0.25 km 2 /s 2 ). 
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Figure 10 . The result of employing simple sail-steering laws for optimally changing the semimajor axis and 
the inclination (C 3 = 0.25 km 2 /s 2 ). The eccentricity is not directly controlled. 
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Figure 11. Three-dimensional trajectory illustration for employing simple sail-steering laws for optimally 
changing the semimajor axis and the inclination. 


available in MATLAB. Downhill-Simplex searches for the initial costates given the boundary conditions. 
The trajectory reaches 0.48 AU with leftover eccentricity, although at some points the eccentricity is very 
small. The period of the eccentricity oscillation is clearly shrinking and the eccentricity is not growing 
but is instead also dropping in magnitude. In both of these ways the result is qualitatively similar to the 
baseline optimal trajectory design by Carl Sauer, indicating that the optimization is capable of controlling 
semimajor axis and eccentricity together but is just not driving the eccentricity down enough. The results are 
a complex combination of the weightings on the final state errors and the final conditions on the numerical 
integration, which currently exits once the target radius is reached. Further work is necessary to determine 
a methodology for selecting the weights and terminating the initial conditions. 

A solar sailing kinetic energy interceptor (KEI) mission, which requires at least ten 160-m, 300-kg solar 
sail spacecraft with a characteristic acceleration of 0.5 mm/s 2 , is proposed in Ref. 12 as a viable near- 
term option for mitigating the threat posed by near-Earth asteroids (NEAs). The solar sailing phase of the 
proposed KEI mission, which is very similar to that of the SPI mission, is comprised of the initial cruise phase 
from 1 AU to 0.25 AU, the cranking orbit phase (for a 168-deg inclination change), and the final retrograde 
orbit phase prior to intercepting the target asteroid at its perihelion. A solar sailing KEI trajectory design 
for optimal intercept /impact /deflection of NEAs, which is applicable to the SPI trajectory design, can be 
found in Ref. 13. 

VI. Attitude Control of the SPI Sailcraft Using TCMs and RSBs 

A. A Simplified Dynamical Model 

Consider a 160-m SPI sailcraft in a heliocentric orbit equipped with trim control masses (TCMs) and roll 
stabilizer bars (RSBs). The roll axis is perpendicular to the sail plane and it often nominally points toward 
the sun (but not necessarily). The pitch and yaw axes are the transverse axes along the masts. A dynamical 
model consists of a rigid sailcraft of mass M and two trim ballasts of each mass of m, as shown in Fig. 13. 
The origin of the body-fixed reference frame is located at point O, which is assumed to be the center of mass 
of a sailcraft of mass M excluding the ballast masses. It is also assumed that TCMs are running along mast 
lanyards (i.e. , along the pitch and yaw axes). 

A set of simplified attitude equations of motion are considered for control design, as follows: 7 


J\U\ - Ui + d\ 

(20) 

J 2^2 — 1 L 2 T d -2 

(21) 

J 3 U 3 = U 3 + d 3 

(22) 
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Figure 12. The result of employing the semimajor axis change optimization using Downhill-Simplex. 
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Figure 13. A solar sail spacecraft with trim control masses in a sun-centered orbit. 
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Figure 14. A dynamic model of sailcraft with a trim control mass moving along the yaw axis to generate a 
pitch control torque. 
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di = ±0.5e.F, di = d 3 = ±eF 


Ji = h + m r (y 2 + z 2 ) 

Ji — h + m r z 2 
Ji = h + m r y 2 

where 5 is the RSB tilt angle, C is the RSB control scale factor, y and 2 are the TCM positions along the 
pitch and yaw axes, respectively, F is the solar radiation pressure force defined as 

F ~ (^ 7 ) F s cos 2 a = 2r] PA cos 2 a (23) 
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and m r is the so-called reduced mass defined as 


m(M + m ) 
M + 2m 


(24) 


m 


Note that m T ~ m because M >> m. 

A reference 160-m SPI sailcraft model for AOCS design is assumed as 

Sail size = 160 x 160 m 
Scallop factor = 75% 

Sail area A = 19,200 m 2 with rj = 0.84 

Solar thrust force F s = 0.16 N (at 1 AU), 0.69 N (at 0.48 AU) 

Sailcraft total mass = 450 kg (150-kg sail, 250-kg bus, and 50-kg payload) 
Characteristic acceleration a c = 0.35 mm/s 2 (at 1 AU) 

Moments of inertia {I\,Iz,h) = (642876, 321490, 321490) kg-m 2 
Cm-cp offset e = 0.4 m (0.25% of 160 m) 

Disturbance torque = 0.064 N-m (1 AU), 0.256 N-m (at 0.48 AU) 

Trim control mass m = 5 kg (each axis) 

Main-body mass M = 440 kg 

TCM speed limit = 5 cm/s 

TCM J/max = Zrnax ~ ±100 m 

RSB maximum deflection angle 5 max = ±45 deg 

RSB moment arm length = 1.7 m 

(a, e, i, = (0.48 AU, IE-5, 45 deg, 0, 0, 0) at t = 0 


B. Roll- Axis Control Design 

The roll-axis control logic is assumed as 


ui = — sat <A'isat[ei + — f e\dt] + D x uji 
Tj J 


(25) 


where u x is the roll control torque command with a saturation limit of t/ max = C sin <5 max and the variable 
limiter L x is self-adjusted as 

L x = ~^ min j i/ 2a il e il> w max } (26) 

where a x = U max / J\ and uj max is assumed to be constrained as 0.1 deg/s. 

The roll control torque u x (in units of N-m) is related to the tilt angle 5 of the RSB of an assumed 
moment arm length of 1.7 m for a 160-m sailcraft, as follows: 


^ . r (1 -7) (53.2) „ . r 
u\ = Csino — Fsmo 


(27) 


It is assumed that all 4 RSBs are rotated simultaneously. Consequently, we have U m &x. of ±0.15 N-m (at 1 
AU) for a maximum tilt angle of ±45 deg. 

The RSB actuator dynamics is also assumed as 


T5 + 5 = <5 C ; |5 c (t)| < <5 n 


(28) 


where T is the actuator time constant, 6 is the actual tilt angle, and 5 C is the commanded tilt angle (the 
control input) with a maximum value of <5 m ax- The commanded tilt angle of the RSB is then given by 


5 C = arcsin | — 


80 u\ 


where 


and u\ is given by Eq. 25. 


7)(53.2)F 

m 


80 u\ 


(1.7)(53.2 )F 


F w 0.16 


cos 2 a 


(29) 

(30) 
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C. Pitch/Yaw Control Design 

The steady-state trim position of the pitch TCM for countering the effect of the cm/cp offset e can be 


estimated as 


M 

z ss = — e ; 
m r 


— e = ±35 m 
m 


The actuator dynamics of the pitch TCM is assumed as 

Ti + 2 = z c ; |z c (t)| < z„ 


(31) 


(32) 


where T is the actuator time constant, z is the actual position, and z c is the commanded position (the control 
input) with a maximum value of Zjnax* for a reference control design, it is assumed that Zm ax — ±100 m, 
i max = ±0.05 m/s, and T - 200 sec for a 160-m sailcraft with M = 440 kg, m = 5 kg, and e = ±0.4 m. 

The pitch control logic is assumed as 


z c = sat 

■Zmax 




4" Z)2^2 


J 


(33) 


where r is the time constant of integral control and variable limiter ±2 is self-adjusted as 


L 2 — jr min|v / 2a 2 |e 2 |, w max | 


Similarly, we have the yaw control logic as 


Vc = ~ sat { K 3 sat[e 3 + — [ e 3 dt\ ± D 3 ui 3 

y max L3 T3 J 


(34) 


(35) 


where y m ax = ±100 m. 

D. Simulation Example: Cranking Orbit Phase at 0.48 AU 

An attitude maneuver for achieving a desired 180-deg <5 clock-angle change within 3 hrs (with a fixed cone- 
angle command of 35 deg) is illustrated in Figs. 15-17. An equivalent, direct 70-deg pitch maneuver within 
3 hrs is also demonstrated in Figs. 18-20. Deviations from the desired cone and clock angles, as can be 
seen in these figures, are caused by the nature of the propellantless attitude control in the presence of 
solar disturbance torques. An inclination increase of 0.06 deg/day can be seen in Fig. 21 for the cranking 
orbit phase of the SPI mission. These preliminary AOCS design results are being validated in using a 
high-fidelity dynamical model of the 160-m sailcraft with actuator nonlinearities and structural flexibility. 11 
The robustness of a propellantless primary ACS described in the preceding section can be further enhanced 
by a secondary ACS utilizing tip- mounted, lightweight pulsed plasma thrusters (PPTs). The microPPT- 
based ACS is mainly intended for attitude recovery maneuvers from off-nominal conditions. Details of the 
microPPT-based ACS design for solar sail spacecraft can be found in Refs. 8 and 14. 

VII. Conclusion 

In this paper, a baseline attitude and orbit control system (AOCS) has been developed for a 160-m, 
450-kg solar sail spacecraft of the Solar Polar Imager (SPI) mission. A reference SPI sailcraft, consisting 
of a 160-m, 150-kg square solar sail, a 250-kg spacecraft bus, and 50-kg science payloads, has a nominal 
solar thrust force of 160 mN (at 1 AU), an uncertain center-of-mass/center-of-pressure offset of ±0.4 m, and 
a characteristic acceleration of 0.35 mm/s 2 . This paper has focused on the solar sailing phase of the SPI 
mission. The simplicity of a quaternion-based attitude control system approach has been demonstrated for 
the case of coupled attitude and trajectory control employing the cone and clock angles of the solar thrust 
vector as trajectory control inputs. A relatively fast, 70-deg pitch reorientation within 3 hrs every half orbit 
during the orbit cranking phase has been shown to be feasible, with the primary ACS, for possible solar 
observations even during the 5-yr cranking orbit phase. 
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Figure 16. A 180-deg S clock-angle maneuver with a fixed 35-deg cone-angle command (continued). 
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Figure 19. A direct 70-deg pitch-axis maneuver (continued). 
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